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Abstract 

We discuss the potential of the Sudbury Neutrino Observatory (SNO) to 
constraint the four-neutrino mixing schemes favoured by the results of all 
neutrino oscillations experiments. These schemes allow simultaneous transi- 
tions of solar u' e s into active z^'s, i/ T 's and sterile v s controlled by the ad- 
ditional parameter cos 2 ($23) cos 2 ($24) and they contain as limiting cases the 
pure i^-active and z^ e -sterile neutrino oscillations. We first obtain the solu- 
tions allowed by the existing data in the framework of the BPOO standard 
solar model and quantify the corresponding predictions for the Charged Cur- 
rent and the Neutral Current to Charged Current event ratios at SNO in the 
different allowed regions as a function of the active-sterile admixture. Our 
results show that some information on the value of cos 2 ($23) cos 2 ("$24) can be 
obtained by the first SNO measurement of the CC ratio, while considerable 
improvement on the knowledge of this mixing will be achievable after the 
measurement of the NC/CC ratio. 
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I. INTRODUCTION 



The Sudbury Neutrino Observatory M is a second generation water Cerenkov detector 
using 1000 tonnes of heavy water, D 2 0, as detection medium. SNO was designed to ad- 
dress the problem of the deficit of solar neutrinos observed previously in the Homestake @ , 
Sage 0, Gallex+GNO [|,|J , Kamiokande || and Super-Kamiokande |7]|| experiments, 
by having sensitivity to all flavours of neutrinos and not just to u e , allowing for a model 
independent test of the oscillation explanation of the observed deficit. 

Such sensitivity can be achievable because energetic neutrinos can interact in the D 2 of 
SNO via three different reactions. Electron neutrinos may interact via the Charged Current 
(CC) reaction 

v e + d — > p + p + e~ , (1) 

with an energy threshold of several MeV. All non-sterile neutrinos may also interact via 
Neutral Current (NC) 

v x + d — > n + p + v' x , (x = e,/i, r), (2) 

with an energy threshold of 2.225 MeV. With smaller cross section the non-sterile neutrinos 
can also interact via Elastic Scattering (ES) v x + e~ — > v' x + e~. 

The main objective of SNO is to measure the ratio of NC/CC events. In its first year 
of operation SNO is concentrating on the measurement of the CC reaction rate while in 
a following phase, after the addition of MgCl2 salt to enhance the NC signal, it will also 
perform a precise measurement of the NC rate. It is clear that a cross-section-normalized 
and acceptance-corrected ratio higher than 1 would strongly indicate the oscillation of v e 
into and/or v T . On the other hand a deficit on both CC and NC leading to a normalized 
NC/CC ratio 1, can only be made compatible with the oscillation hypothesis if v e oscillates 
in to a sterile neutrino. 

There are several detailed studies in the literature of the potential of the SNO experiment 
to discriminate between the different oscillation solutions to the solar neutrino problem 
(SNP) P|-|r2j|. Most of these studies have been performed in the framework of oscillations 
between two neutrino states where v e oscillates into either an active, v e — > u a , or a sterile, 
v e ~~ > v si neutrino channel. On the other hand, once the possibility of a sterile neutrino is 
considered, these two scenarios are only limiting cases of the most general mixing structure 
13| , |14|1 which permits simultaneous v e — > v s and v e — > v a oscillations. 



In this paper we study the potential of the Sudbury Neutrino Observatory to discriminate 
between active and sterile solar neutrino oscillations when analyzed in the framework of four- 
neutrino mixing. We consider those four-neutrino schemes favoured by considering together 
with the solar neutrino data, the results of the two additional evidences pointing out towards 
the existence of neutrino masses and mixing: the atmospheric neutrino data [JTSJ and the 
LSND results [pl|. We concentrate on two SNO measurements: the first expected result on 



the CC ratio and the expected to be most sensitive, the ratio of NC/CC. The measurement 
of other observables, such as the recoil energy spectrum of the CC events and the zenith 
angular dependence |9Hl2| can provide important information to distinguish between the 
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different allowed regions for i/ e -active oscillations but they are not expected to be very 
sensitive as discriminatory between the active and sterile oscillations. 

The outline of the paper is the following. For the sake of completeness we begin by 
discussing in Sec. |I| the expected results when obtained in the pure two-neutrino oscillation 
hypothesis. In Section |TTT| we determine the presently allowed regions for the oscillation 
solutions to the SNP in the framework of four-neutrino mixing. In Sec. |IV| we present 
the results of the expected CC and NC/CC rates for the different solutions and quantify 
the attainable sensitivity to the additional mixing controlling the admixture of active-sterile 
entering into the solar neutrino oscillations. Finally in Sec. |V| we summarize our conclusions. 



II. TWO NEUTRINO MIXING: ALLOWED REGIONS AND PREDICTIONS 

FOR SNO 

We first describe the results of the analysis of the solar neutrino data in terms of v e 
oscillations into either active or sterile neutrinos. We determine the allowed range of oscil- 
lation parameters using the total event rates of the Chlorine @, Gallium and Super- 
Kamiokande (corresponding to the 1117 days data sample) experiments. For the Gal- 
lium experiments we have used the weighted average of the results from GALLEX+GNO 
and SAGE detectors. We have also include the Super-Kamiokande electron recoil energy 
spectrum measured separately during the day and night periods. This will be referred in 
the following as the day-night spectra data which contains 18 + 18 data bins. The analysis 



includes the latest standard solar model fluxes, BPOO model fll?! , with updated distributions 
for neutrino production points and solar matter density. For details on the statistical analy- 
sis applied to the different observables we refer to Ref. |T8| , [T9|j . Nevertheless, two comments 
on the statistical analysis are in order: 

- In the present analysis we also include the contribution to the theoretical errors of the 
event rates arising from the small uncertainty in the measured So-factor for the reaction 
16 0(p,7) 17 F which is new in the BPOO model as discussed in Ref. JTTj] . Following the standard 
procedure p0| , we include this new source of uncertainty for the rates, that we denote as Cp, 
by adding a new fractional la uncertainty AlnXc F = 0.18. Since this uncertainty affects in 
direct proportion to the 17 F flux we correspondingly add a new line a>ic F = d§i/dX CF to 
the response matrix, with values a ^f,c f = 1 an d ctic F — for all other fluxes. 
- In the analysis of the day-night spectrum data we include the correlation between the 



systematic errors of the day and night bins which were conservatively ignored in Ref. |L9 
Thus, we use the correlation matrix: 



a 



ij ~ ^ij( a i,stat a i,uncorr) + a i,exp& j,exp + Vi^cal&j^al (3) 



where i and j run from 1 to 36 bins in the day-night spectra data. a itSt at is the statistical 
error, and a ijUncorr is the error due to uncorrelated systematic uncertainties. a i)exp and a itCa i 
are the correlated errors due to correlated systematic experimental uncertainties and the 
calculation of the expected spectrum respectively (see Ref. |L8| for details). The addition of 



the correlations between the errors for the day and night bins, which more properly takes 
into account the day-night information, leads to stronger constraints on the regeneration 
region. 
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With all this we obtain that using the predicted fluxes from the BPOO model the \ 2 f° r 
the total event rates is Xssm = 56 for 3 d. o. f. This means that the SSM together with 
the SM of particle interactions can explain the observed data with a probability lower than 
5 x 1(T 12 . 

The allowed regions in the oscillation parameter space are shown in Fig. [I]. We present 
them in the full parameter space for oscillations including both MSW |21| and vacuum | 22| 



oscillations, as well as quasi-vacuum |E3| oscillations (QVO) and matter effects for mixing 



angles in the second octant (the so called dark side ||^ , |T4l , |T9|| ) . In the case of i/ e -active 
neutrino oscillations we find that the best-fit point is obtained for the LMA solution. There 
are two more local minima of \ 2 in the MSW region: the SMA and LOW solutions. Notice 
also that LOW and QVO regions are connected at the 99 %CL and they extend into the 
second octant so maximal mixing is allowed at 99 % CL for Am 2 in what we define as the 
LOW-QVO region. 

Following the standard procedure, the allowed regions are defined in terms of shifts of 
the x 2 function with respect to the global minimum in the plane. Defined this way, the size 
of a region depends on the relative quality of its local minimum with respect to the global 
minimum but from the size of the region we cannot infer the actual absolute quality of the 
description in each region. In order to give this information we list in Table | the goodness 
of the fit (GOF) for each solution obtained from the value of \ 2 at the different minima. 

For oscillations into sterile neutrinos the global minimum lies in the SMA solution. As 
seen in Fig. p] we find that with the present data and using the criteria explained above, 
there are also allowed solutions for sterile neutrinos in the LMA and LOW-QVO regions 
at 99 %CL once the day-night spectra data is included. We consider, however, that they 
are not acceptable solutions as their fit to the global rates is really poor with a probability 
of acceptance less than 0.004 []. The differences between both oscillation scenarios (active 
and sterile) can be easily understood. Unlike active neutrinos which lead to events in the 
Super-Kamiokande detector by interacting via NC with the electrons, sterile neutrinos do 
not contribute to the Super-Kamiokande event rates. Therefore a larger survival probability 
for 8 B neutrinos is needed to accommodate the measured rate. As a consequence a larger 
contribution from 8 B neutrinos to the Chlorine and Gallium experiments is expected, so 
that the small measured rate in Chlorine can only be accommodated if no 7 Be neutrinos 
are present in the flux. This is only possible in the SMA solution region, since in the LMA 
and LOW regions the suppression of 7 Be neutrinos is not enough. Notice also that the SMA 
region for oscillations into sterile neutrinos is slightly shifted downwards as compared with 
the active case. This is due to the small modification on the neutrino survival probability 
induced by the different matter potentials. The matter potential for sterile neutrinos is 
smaller than for active neutrinos due to the negative NC contribution proportional to the 
neutron abundance. For this reason the resonant condition for sterile neutrinos is achieved 
at lower Am 2 . On the other hand, the flatter spectrum, is better fitted in both LMA and 
LOW regions independently of the active or sterile nature of the neutrino. This leads to the 



* Marginally allowed VO solutions were also possible (see for instance Ref. ||25| ) with last year data 
sample but they are now ruled out. 
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improvement of the quality of the description for these solutions for both active and sterile 
neutrinos. However, as mentioned above, for the analysis of the total rates these LMA and 
LOW solutions give a very bad fit in the sterile case and we decide not to consider them 
in the following. Also, as we will see in next section, when the analysis is performed in the 
framework of four-neutrino oscillations those large mixing solutions for sterile neutrinos do 
not appear. 

Next we quantify the predictions for the SNO observables in the allowed regions discussed 
above. The total number of events in the CC reaction at SNO can be obtained as 



N&c = E « 

fe=l,2 



JdE u \ k (E u )a cc (E u )(P u ^ Ue ) (4) 



where E v is the neutrino energy, <pk are the total neutrino 8 B and hep fluxes, is the neutrino 
energy spectrum (normalized to 1) and {P Ve -*v e ) is t ne time-averaged v e survival probability 
for oscillations into either active or sterile neutrinos. Here occ is the vd CC cross section 
computed from the corresponding differential cross sections folded with the finite energy 
resolution function of the detector and integrated over the electron recoil energy: 

°oo{Eu) = j T JT J dT> Res(T, T') dffcg ffi T) , (5) 

where T and T' are the measured and the true kinetic energy of the recoil electrons and T t h 
indicates the threshold expected from the experiment. The resolution function Res(T, T') 
is of the form ||: 

1 " (T — T^ 2 

Res(T, T ) = - exp 

v / 2^(0.348jT'/MeV) 



0.242 T> MeV 



(6) 

and we take the differential cross section dacc(E u , T')/dT' from p6| . For definiteness,we 
adopted the most optimistic total energy threshold E th = 5 MeV (T th = E th — m e ). 

Correspondingly, the total number of events in the NC reaction at SNO is obtained as 

Kc = E fa j dE » h(E v )<T NC (E u )((P Ve ^ e ) + (P„ e ^ » (7) 

k=l,2 J 



where a^c is the vd NC cross section from |26| and (P Ue ^ Ua ) is the time-averaged probability 



of oscillation into any other active neutrino. In the case that v e oscillates only into active 
neutrinos (P Ve -*v e ) + (Pu e ^u a ) = ^ and N$ c is a constant. 

In order to cancel out all energy independent efficiencies and normalizations we will use 
the ratio: 

th Ncc 

R CC = N SSM = t CC ] ( 8 ) 

where N^q M is the predicted number of events in the case of no oscillations. The equivalent 
expression for the NC ratio 

j\rth 

Rnc = = [ NC ] ( 9 ) 
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Out of those ratios one can compute the double ratio -$p- = NC / [CC for which the largest 

H cc 



sources of uncertainties cancel out |TT| . As it was shown in Ref . |26[ , the ratio between the 
NC and CC reaction cross sections is extremely stable against any variations of the inputs 
of the calculations. The expected total uncertainties for the [CC] ratio and the [NC]/[CC] 
ratio are 6.7 % and 3.6 % respectively assuming 5000 CC events and 1219 NC events |Ljfl . 

In Fig. we show the predicted [CC] and [NC]/[CC] ratios for the allowed regions in 
the two flavour analysis. The dots correspond to the local best fit points and the error bars 
show the range of predictions for the points inside the 90 and 99 %CL allowed regions. The 
mapping of the regions onto these bars can be easily understood from the behaviour of the 
probability for the different solutions: 

(a) For oscillations into active neutrinos the [NC]/[CC] ratio is simply the inverse of the 
[CC] prediction. Therefore 

• In the SMA region smaller mixing angles are mapped onto higher (lower) values of 
[CC] ([NC]/[CC]) ratio. One may notice that the prediction for the [CC] rate for the 
global best fit point (0.72) is larger than the measured rate at Super-Kamiokande. 
This is due to the nearly flat spectrum at Super-Kamiokande which implies that the 
best fit point in the global analysis corresponds to a smaller mixing angle than the 
best fit point for the analysis of rates only. 

• In the LMA region, the lower Am 2 and 9 values are mapped onto higher (lower) 
[NC]/[CC] ([CC]) ratios and viceversa. 

• In the LOW region the higher (lower) [NC]/[CC] ([CC]) ratio occurs for smaller 9 and 
higher Am 2 . 

(b) For the sterile case, the best fit point in SMA occurs at lower Am 2 than in the active 
case and this produces a higher prediction for the [CC] ratio (0.74). The [NC]/[CC] ratio 
takes an almost constant value very close to one (0.98 in the best fit point), since both 
numerator and denominator are proportional to (P Ue ^ Ue ). It is smaller than one because for 
the SMA solution the probability increases with energy in the range of detection at SNO 
and the threshold for the NC reaction is below the one for the CC one. 

For the sake of consistency we have checked that our results agree perfectly with those in 



Ref. ]TT|] when comparing the same points in the parameter space. However a careful reader 
may notice that the predictions at the best fit points and ranges in each region displayed in 
Fig. are slightly different of those in Ref. [pjj . The difference is due to two factors. First, 



the allowed regions are defined in a different way. In Ref. |TTJ departures from the standard 



solar model in the boron flux normalization are allowed and moreover the regions are defined 
in terms of shifts of the \ 2 function with respect to the local minimum in the corresponding 
region. Second, the inclusion of the updated data, mainly the Super-Kamiokande day-night 
spectra, lowers the value of tan 2 9 for the best bit point in the SMA region by a factor 2 
and increases Am 2 for the best fit point in LMA by a factor 1.5. 

What we see from these results is that while the data on [CC] can give a hint towards 
large or small mixing solutions, it will be hard to distinguish active from sterile oscillations 
on the only bases of this measurement. This is not the case for the [NC]/[CC] ratio where 
both scenarios appear nicely separated. It is not hard to foresee from these results that from 
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the [NC]/[CC] measurement SNO will be able to constraint the additional mixings in the 
four-neutrino scenario which describe the admixture of active and sterile oscillations. This 
is the main point in this paper. 



III. ALLOWED FOUR NEUTRINO MIXING PARAMETERS 

Together with the results from the solar neutrino experiments we have two more evidences 
pointing out towards the existence of neutrino masses and mixing: the atmospheric neutrino 



data |15j and the LSND results [jig] . All these experimental results can be accommodated 
in a single neutrino oscillation framework only if there are at least three different scales of 
neutrino mass-squared differences. The simplest case of three independent mass-squared 
differences requires the existence of a light sterile neutrino, i.e. one whose interaction with 
standard model particles is much weaker than the SM weak interaction, so it does not affect 
the invisible Z decay width, precisely measured at LEP. 

There are six possible four-neutrino schemes that can accomodate all these evidences. 
They can be divided in two classes: 3+1 and 2+2. In the 3+1 schemes there is a group of 
three neutrino masses separated from an isolated mass by a gap of the order of leV which 
gives the mass-squared difference responsible for the short-baseline oscillations observed 
in the LSND experiment. In 2+2 schemes there are two pairs of close masses separated 
by the LSND gap. We have ordered the masses in such a way that in all these schemes 
Am s 2 un = Amji produces solar neutrino oscillations and AmL SND = Am| x (we use the 
common notation Am^ = mf. — m|). 3+1 schemes are disfavoured by experimental data 
with respect to the 2+2 schemes P7|,PS| but they are still marginally allowed [E5|. 



In any of these four-neutrino schemes the flavour neutrino fields v a i (we choose a = 
e, s, n, t) are related to the fields u^l of neutrinos with masses by a rotation U. U is a 4x4 
unitary mixing matrix, which contains, in general, 6 mixing angles and 3 CP violating phases 
(3 additional phases appear for Majorana neutrinos but they are irrelevant for oscillations). 
We neglect here the CP phases, which, in the schemes considered, are irrelevant for solar 
neutrinos because their effect is washed out by averaging over neutrino energy and distance. 
Existing bounds from negative searches for neutrino oscillations performed at colliders as 
well as reactor experiments, in particular the negative results of the Bugey |30| and CHOOZ 



31j v e disappearance experiment, impose severe constrains on the possible mixing structures 
for the four-neutrino scenario. In particular they imply that the matrix elements and 
U e 4 are very small p7| , p8| , |32l . As a consequence, for any of these four-neutrino schemes, 
either 2+2 or 3+1, only four mixing angles are relevant in the study of solar neutrino 
oscillations |l3| , p!4]j32]| and the U matrix can be written as 



U 



( Cia S12 ^ 

— S 12C23C24 C12C23C24 S23C24 S24 

■S12(C23"S24"S34 + S23C34) — Ci2(.S23 c 34 + C23S24S34) c 23 c 34 — S23S24S34 C 24-S34 

V ■Si2(c23S24C34 — S23S34) C\2 (S23 s 34 — C23S24C34) — (C23S34 + S23S24C34) C24C34 / 



(10) 



where $12, $23, $24, $34 are four mixing angles and Cy = cos^jj and Sij = sint^-. 

Since solar neutrino oscillations are generated by the mass-square difference between i>2 
and v\ , it is clear from Eq. (|Il]) that the survival of solar v e 's mainly depends on the mixing 
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angle $12, whereas the mixing angles $23 and $24 determine the relative amount of transitions 
into sterile v s or active v a , this last one being a combination of and v T controlled by the 
mixing angle #34. and v T cannot be distinguished in solar neutrino experiments, because 
their matter potential and their interaction in the detectors are equal, due only to NC weak 
interactions. As a consequence the active/sterile ratio and the survival probability for solar 
neutrino oscillations do not depend on the mixing angle $34, and depend on the mixing angles 
$23 $24 only through the combination cos $23 cos $24- For further details see Ref. [I3"|,ri|. We 
distinguish the following limiting cases: 

• cos $23 cos $24 = corresponding to the limit of pure two-generation v e — > v a transitions. 

• cos $23 cos $24 = 1 for which we have the limit of pure two-generation v e — > v s transitions. 

• If cos $23 cos $24 7^ 1, solar UgS can transform in the linear combination v a of active and 
v T . 

In the general case of simultaneous v e — > v s and v e — > v a oscillations the corresponding 
probabilities are given by [0.0] 

Pu e ^u s — C 23 C 24 (1 _ Pv e -»/e) J (11) 
P^a = (l - C1 3 C 2 24 ) (1 - P^J . (12) 

where P Ve -t Ve takes the standard two-neutrino oscillation form for Am 2 2 and #12 but com- 
puted with the modified matter potential 

A = A cc + cj 3 c 2 M A NC . (13) 

Thus the analysis of the solar neutrino data in the four-neutrino mixing schemes is equivalent 
to the two-neutrino analysis but taking into account that the parameter space is now three- 
dimensional (Am^ 2 , tan 2 1? 12 , cos 2 $23 cos 2 $24). We want to stress that, although originally 



this derivation was performed in the framework of the 2+2 schemes [pT3l,p~4| , it is equally 
valid for the 3+1 ones [32| . 



We first present the results of the allowed regions in the three-parameter space for the 
global combination of observables. Notice that since the parameter space is 3-dimensional 
the allowed regions for a given CL are defined as the set of points satisfying the condition 
X 2 (Am 2 2 ,tf 12 ,c 2 3C 2 4) - xLn < AX 2 (CL, 3 dof) where, for instance, A X 2 (CL, 3 dof)=6.25, 
7.83, and 11.36 for CL=90, 95, and 99 % respectively. In Figs. |]we plot the sections of such 
volume in the plane (Am^, tan 2 (-$ 12 )) for different values of C2 3 c 2 4 . The global minimum used 
in the construction of the regions lies in the LMA region and for pure i/ e -active oscillations, 

c 23 c 24 — u - 

As seen in Fig. [3] the SMA region is always a valid solution for any value of c 2 3 c 2 4 at 99% 
CL (the same is true at 95% CL). This is expected as in the two-neutrino oscillation picture 
this solution holds both for pure ^ e -active and pure z/ e -sterile oscillations. Notice, however, 
that the statistical analysis is different: in the two-neutrino picture the pure z/ e -active and 
z/ e -sterile cases are analyzed separately, whereas in the four-neutrino picture they are taken 
into account simultaneously in a consistent scheme. Since the GOF of the SMA solution 
for pure i/ e -sterile oscillations is worse than for SMA pure active oscillations (as discussed 
in the previous section), the corresponding allowed region is smaller because they are now 
defined with respect to a common minimum. Also, we notice, that for the SMA solution the 
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best scenario is a non-zero admixture between active and sterile oscillations. For this reason 
this solution is allowed at a CL better than 90 % only in the range 0.11 < c| 3 c| 4 < 0.31. 

On the other hand, the LMA and LOW-QVO solutions disappear for increasing values 
of the mixing c^c^. We list in Table |T|the ranges of c^%c\ A for which each of the solutions 
is allowed at a given CL. We see that at 95 %CL the LMA solution is allowed for maximal 
active-sterile mixing c\%c\ A = 0.5 while at 99%CL all solutions are possible for maximal 
admixture. 



IV. EXPECTED RATES AT SNO IN FOUR NEUTRINO SCHEMES 

In this section, we present the predictions for the CC ratio and for the NC/CC ratio in 
the four-neutrino scenario previously described. This scenario contains as limiting cases the 
pure z/ e -active and z/ e -sterile neutrino oscillations. However, when comparing the results 
for both limiting cases with the ones presented in Sec. |TJ the reader must notice that there 
are some changes in the predicted ranges because the allowed regions are obtained with a 
different statistical criteria. Now, as discussed above, all the allowed regions are defined 
with respect to the same global minimum (laying in the LMA with c\^c\ A =Q) with 3 dof. 
Because of that, the predicted ranges in the four-neutrino scheme are wider for the pure 
i/ e -active oscillations and narrower for the v e -sterile case. 

In Figs. £|— 16| we show the results for the predicted [CC] ratio and [NC]/[CC] ratio for 
the different allowed regions (SMA, LMA, LOW-QVO) at 90 and 99 %CL as a function of 
c 23 c 24- The general behaviour of the dependence of the predicted ratios with c^c\ A can be 
easily understood using the following simplified expressions obtained from Eqs. @, @ and 

(0): 

[CC] ~ P n , (14) 
[NC] 1 - c| 3 cl 4 (l - P^J 



[CC] p 



(15) 



V e - 



From Eq. (|H|) we see that the only dependence of [CC] on c| 3 c| 4 is due to the modification of 
the matter potential entering in the evolution equation (see Eq. (|H|) and discussion above) 
and it is very weak. The dependence of the allowed range of the [CC] ratio with c^c^ 
displayed in the figures arises mainly from the variation of the size of the allowed regions. 
Alternatively following Eq. ( |T5D we find a stronger linear dependence of [NC]/[CC] on c^c^ 
with slope (-1+P^J/P^ ~ 1-1/[CC] and intercept l/P Ve -> Ve ~ 1/[CC]. This simple 
description is able to reproduce the main features of our numerical calculations as can be 
seen in the figures. 

Figs. f|— P contain the main quantitative result of our analysis in the four-neutrino mix- 
ing scenario. From each of them it is possible to infer the allowed range of the active-sterile 
admixture, c^c^, compatible, within the expected uncertainty, with a given SNO mea- 
surement of the ratios. Also, comparing the allowed ranges for the different solutions one 
can study the potential of these measurements as discriminatory among the three presently 
allowed regions. Of course, both issues are not independent as we have no "a priory" knowl- 
edge of which is the right solution and both must be discussed simultaneously. In order to 
do so we pass to describe and compare in detail the predictions in the different regions. 
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The results for the SMA solution are shown in Fig. |.a and g.b for [CC] and [NC]/[CC] 
ratios respectively. First we notice that we find a small region allowed at 90 %CL only for a 
non-vanishing admixture of active and sterile oscillations as mentioned before. In this region 
[CC]~0.65-0.73 and [NC]/[CC]~ 1.3-1.4. The predictions at 99% range from [CC]~ 0.4-0.9 
([NC]/[CC]~1. 1-2.5) for pure z/ e -active scenario to to [CC]~ 0.59-0.85 ([NC]/[CC]~0.96- 
0.98) for pure z/ e -sterile oscillations. Thus if SNO observes a ratio [CC]<0.58 the value of 
c 23 c 24 can be constrained to be smaller than 1 disfavouring pure z/ e -sterile oscillations. On 
the contrary a measurement of [CC]>0.68 will immediately hint towards the SMA solution 
but will not provide any information on the active-sterile admixture. Also, one must notice, 
that such value, although allowed by the present global statistical analysis at 99 %CL, will 
imply a strong disagreement with the total rate event rate observed at Super-Kamiokande. 

As seen in Fig. f|.b the [NC]/[CC] ratio is more sensitive to the active-sterile admixture. 
To guide the eye, in the figures for the [NC]/[CC] ratio we plot a dotted line for the prediction 
in the case of no oscillation [NC]/[CC]=1. For any of the solutions, the allowed range for 
this ratio shows as general behaviour a decreasing with c| 3 c| 4 due to two effects: (i) the 
allowed regions become smaller and (ii) the prediction decreases when more sterile neutrino is 
involved in the oscillations as described in Eq. (|15j). The measurement of higher values of this 
ratio will favour the four-neutrino scenario with larger component of z/ e -active oscillations. 
On the other hand a measurement of [NC]/[CC]~ 1, will push the oscillation hypothesis 
towards the pure z/ e -sterile oscillation scenario. This case will be harder to differentiate 
from the non-oscillation scenario. We find that with the expected sensitivity the parameter 
c 23 c 24 i s constrained to be above 0.44 at 99% CL and that the pure z/ e -active oscillations in 
the SMA region are compatible with [NC]/[CC]= 1 only at ~ 5a. 

The predictions for oscillation parameters in the LMA region are shown in Figs. |5|.a 
and |5].b for [CC] and [NC]/[CC] ratios respectively. The predictions at 99% vary in the 
range [CC]~ 0.18-0.62 and [NC]/[CC]~ 1.4-5.6. The first thing we notice by comparing 
Fig. [|.a with Fig. ^.a and Fig. |^.a is that the most discriminatory scenario for the [CC] rate 
results if SNO finds a small value [CC]~ 0.25. This would significantly hint towards the 
LMA solution to the solar neutrino problem and towards and z/ e -active oscillation scenario. 
First, it is well separated from the predictions for the SMA and LOW regions. Second, it will 
include as a bonus a small but measurable day-night asymmetry [|H],|ll]]. And third it will 
constrain the C2 3 c| 4 to a small value (~ 0.2). On the contrary the less discriminatory scenario 
will be a measurement 0.4<[CC]<0.6 where the prediction would be compatible with both 
SMA and LOW-QVO solutions and no improvement on our knowledge of the four-neutrino 
schemes is possible. The [NC]/[CC] ratio can definitively improve the discrimination between 
the different scenarios provided its measurement lies in the upper range. For instance a 
measurement of [NC]/[CC]~ 4 (± 0.7 at 5a) will be conclusive for selecting LMA as the 
solution to the SNP and will imply and upper bound on ci^c^ < 0.3. 

The predictions for the LOW-QVO region lie between the ones for SMA and LMA as 
displayed in Fig. || and therefore they are more difficult to discriminate. The predictions at 
99% vary in the range [CC]~ 0.3-0.68 and [NC]/[CC]~ 1.2-3.4. As a consequence we see that 
a low [CC] ratio but still within the 99% CL range allowed for this region, 0.3<[CC]<0.4, 
will constrain significantly the c^c^ parameter compatible with this solution but it will not 
be distinguishable from the LMA solution unless the measured [CC]<0.3. As mentioned 
above, the [NC]/[CC] ratio will be able to differentiate the LMA and LOW-QVO solutions 
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if not in the range [1.5,3]. One should also notice that for the upper part of this range a 
positive measurement of the day-night asymmetry and the zenith dependence [12| will point 
towards the higher Am 2 of the LOW region as the solution. 



V. DISCUSSION 

In this paper we have studied the potential of the Sudbury Neutrino Observatory to dis- 
criminate between active or sterile solar neutrino oscillations when analyzed in the framework 
of four-neutrino mixing. We considered those four-neutrino schemes favoured by consider- 
ing together with the solar neutrino data, the results of the two additional evidences pointing 



out towards the existence of neutrino masses and mixing: the atmospheric neutrino data [15 



and the LSND results (TIJ. These schemes allow simultaneous transitions of solar v' e s into 



active z/^'s, z/ T 's and sterile v s controlled by the additional parameter cos ($23) cos ($24) and 
they contain as limiting cases the pure z/ e -active and z/ e -sterile neutrino oscillations. The 
allowed solar solutions have been reanalyzed including the recently BPOO standard solar 
model and the latest solar neutrino data. We find that the global minimum lies in the LMA 
region and for pure t/ e -active oscillations {c\^c\ A = 0). We also find that in the framework of 
four-neutrino mixing the SMA solution is allowed at 90% CL for non vanishing active-sterile 
mixing c^c 2 ^ in the range [0.11,0.31]. 

We concentrated on two SNO measurements: the first expected result on the [CC] ratio 
and the expected to be most sensitive to the active-sterile admixture, the ratio of [NC]/[CC] 
and evaluated the predictions in the different regions as a function of the additional mixing 
c 23 c 24- Our results are display in Figs. [|— [| They show that in most cases with the mea- 
surement of the [CC] ratio, it will be hard to improve the present knowledge of cf^c 2 ^ but 
with the precise determination of the [NC/CC] ratio at SNO, this parameter can be strongly 
constrained for some of the allowed solutions. For example, we find that for the [CC] rate 
the most discriminatory scenario would be that SNO finds a small value [CC]~ 0.25. This 
significantly hints towards the LMA solution to the solar neutrino problem and towards 
an active-active oscillation scenario. In this case the [NC]/[CC]~ 4 (± 0.7 at 5cr) will be 
conclusive for selecting LMA as the solution to the SNP and will imply and upper bound on 
c 23 c 24 < 0-3- Conversely, a measurement of [NC]/[CC]~ 1, although harder to distinguish 
from the non-oscillation scenario, will push the oscillation hypothesis towards the sterile 
SMA solution, and with the expected sensitivity a bound c^c^ > 0.44 at 99% CL can be 
imposed. 



ACKNOWLEDGMENTS 

We thank J. N. Bahcall for clarifying several aspects of BPOO standard solar model. 
We also acknowledge valuable discussions with E. Lisi on the correlations in the day-night 
spectra and with C. Giunti on the four-neutrino schemes. This work was supported by the 
Spanish DGICYT under grants PB98-0693 and PB97-1261, by the Generalitat Valenciana 
under grant GV99-3-1-01 and by the TMR network grant ERBFMRXCT960090 of the 
European Union. 



10 



REFERENCES 



[1] A. B. McDonald, Nucl. Phys. B (Proc. Suppl.) 77, 43 (1999); SNO Collaboration, 
Physics in Canada 48, 112 (1992); SNO Collaboration, |nucl-ex/99100T6 . 

[2] B. T. Cleveland et al, Astrophys. J. 496, 505 (1998); R. Davis, Prog. Part. Nucl. Phys. 
32, 13 (1994). 

[3] SAGE Collaboration, J. N. Abdurashitov et al, Phys. Rev. C60, 055801 (1999); Talk 
by V. Gavrin at Neutrino 2000, Sudbury, Canada, June 2000 Qhttp: / /nu2000.sn~ 
laurentian.ca|) . 



[4] GALLEX Collaboration, W. Hampel et al, Phys. Lett. B447, 127 (1999). 
[5] Talk by E. Belloti at Neutrino 2000, Sudbury, Canada, June 2000 ( |http://nu2000.sno.- 
laurentian.cap . 

[6] Kamiokande Collaboration, Y. Fukuda et al, Phys. Rev. Lett. 77, 1683 (1996). 

[7] Super-Kamiokande Collaboration, Y. Fukuda et al, Phys. Rev. Lett. 81, 1158 (1998); 

Erratum 81, 4279 (1998); 82, 1810 (1999); 82, 2430 (1999); Y. Suzuki, Nucl. Phys. B 

(Proc. Suppl.) 77, 35 (1999). 
[8] Talk by Y. Suzuki at Neutrino 2000, Sudbury, Canada, June 2000 ( [http: / / nu2000.sn~ 

laurentian.ca|) . 



[9] J. N. Bahcall and E. Lisi, Phys. Rev. D54, 5417 (1996). 
[10] G. L. Fogli, E. Lisi, D. Montanino, Phys. Lett. B 434, 333 (1998); W. Kwong and 
S. P. Rosen, Phys. Rev. D 54 2043 (1996); F.L. Villante, G. Fiorentini, and E. Lisi, 
Phys. Rev. D 59 013006 (1999);G. L. Fogli, E. Lisi, D. Montanino, and A. Palazzo, 
Hiep-ph/00080l3 ; J- N. Bahcall, P. I. Krastev, A. Yu. Smirnov, Phys. Lett. B 477, 401 
(2000); M. Maris and S. T. Petcov, Phys. Rev. D 62 93006 (2000). 
[11] J. N. Bahcall, P. I. Krastev, A. Yu. Smirnov, Phys. Rev. D62 093004 (2000); [hep 

| ph/0006078| . 

[12] M. C. Gonzalez-Garcia, C. Pefia-Garay, Y. Nir and A. Yu Smirnov, |hep-ph/ 0007227 . 

M. C. Gonzalez-Garcia, C. Pena-Garay and A. Yu Smirnov, in preparation. 
[13] D. Dooling, C. Giunti, K. Kang and C.W. Kim, Phys. Rev. D61, 073011 (2000). 
[14] C. Giunti, M. C. Gonzalez-Garcia, C. Pefia-Garay, Phys. Rev. D62, 013005 (2000). 
[15] NUSEX Collab., M. Aglietta et al., Europhys. Lett. 8, 611 (1989); Frejus Collab., Ch 
Berger et al, Phys. Lett. B227, 489 (1989); 1MB Collab., D. Casper et al, Phys 
Rev. Lett. 66, 2561 (1991); R. Becker-Szendy et al., Phys. Rev. D46, 3720 (1992) 
Kamiokande Collab., H. S. Hirata et al., Phys. Lett. B205, 416 (1988) and Phys. Lett 
B280, 146 (1992); Kamiokande Collab., Y. Fukuda et al, Phys. Lett. B335, 237 (1994) 
Soudan Collab., W. W. M Allison et al, Phys. Lett. B391, 491 (1997). 
[16] C. Athanassopoulos, Phys. Rev. Lett. 75 2650 (1995); Phys. Rev. Lett. 77 3082 (1996); 

Phys. Rev. Lett. 81 1774 (1998). 

[17] J. N. Bahcall, S. Basu and M. H. Pinsonneault, |astro-ph /00 1 0346] , See also, 

[http:/ /www.sns.ias.edu/^bahcall 
[18] M. C. Gonzalez-Garcia, P. C. de Holanda, C. Pefia-Garay, J. W. F. VaUe, Nucl. Phys. 

B573, 3 (2000). 
[19] M. C. Gonzalez-Garcia, C. Pefia-Garay, |hep-ph/000904l| . 

G. L. Fogli, E. Lisi, D. Montanino, and A. Palazzo, Phys. Rev. D62, 013002 (2000). 



11 



[21] S.P. Mikheyev and A.Yu. Smirnov, Sov. Jour. Nucl. Phys. 42, 913 (1985); L. Wolfenstein, 

Phys. Rev. D17, 2369 (1978). 
[22] V.N. Gribov and B.M. Pontecorvo, Phys. Lett. 28B, 493 (1969); V. Barger, K. Whis- 
nant, R.J.N. Phillips, Phys. Rev. D24, 538 (1981); S.L. Glashow and L.M. Krauss, 
Phys. Lett. 190B, 199 (1987). 
[23] S. Petcov Phys. Lett. B214 139(1988) Phys. Lett. B224 426 (1989); J. Pantaleone, 
Phys. Lett. B251,618 (1990); S. Pakvasa and J. Pantaleone, Phys. Rev. Lett. 65, 2479 
(1990); A. Friedland, Phys. Rev. Lett. 85, 936 (2000). 
[24] G. L. Fogli, E. Lisi and D. Montanino, Phys. Rev. D54, 2048 (1996); A. de Goueva, 
H. Murayama and A, Friedland, Phys. Lett. B490, 125 (2000); M. C. Gonzalez-Garcia, 
C. Peha-Garay, Phys. Rev. D62, 031301 (2000). 
[25] V. Barger and K. Whisnant, Phys. Lett. B456 54 (1999); S. Goswani, D. Majumdar 

and A. Raychaudhury, [hep-ph/9909453 . 
[26] S. Nakamura, T. Sato, V. Gudkov and K. Kubodera, frrucl-th/ 00090 12| . See also, 

|http:/ /nuc003.psc.sc.edu/^kubodera 
[27] S.M. Bilenky, C. Giunti and W. Grimus, Eur. Phys. J. C 1, 247 (1998), Proc. of Neutrino 
'96, Helsinki, June 1996, edited by K. Enqvist et ai, p. 174, World Scientific, 1997, |hcp- 
ph/ 9609343| ; S. M. Bilenky, C. Giunti, W. Grimus, and T. Schwetz, Phys. Rev. D60, 



073007 (1999). 

[28] V. Barger, S. Pakvasa, T. J. Weiler, and K. Whisnant, Phys. Rev. D58, 093016 (1998). 

J. J Gomez-Cadenas and M. C. Gonzalez-Garcia, Z. Phys. C71, 443 (1996). 
[29] V. Barger, B. Kayser, J. Learned, T. Weiler, K. Whisnant, Phys. Lett. B489, 345 

(2000); O. L. G. Peres , A. Yu. Smirnov, |hep-ph/0011054 
[30] B. Achkar et ai., Nucl. Phys. B424, 503 (1995). 
[31] M. Apollonio, et ai., Phys. Lett. B466, 415 (1999). 
[32] C. Giunti and M. Laveder, |hep-p.j/0010009| . 



12 



FIGURES 



M ■ 

£ 



io J 



to~ 



10' 



to- 



io-' 



BPOO 




Active 
Global =Rates+Sp[,+|Sp N 

_i I i I i I i I i 



"i 1 ■ i rr 



BPOO 



Sterile 

Global = Rates +Sp p +jSp M 



i i i i i 



IV* 



10' 



iv 1 to- 1 t 



10~ 



to' 



1 to. 
tan 2 i? 



FIG. 1. 90, 95 and 99 %CL allowed regions in the two-neutrino oscillation scenario from the 
global analysis of solar neutrino data including the total measured rates and the Super-Kamiokande 
measured spectrum at day and night. The global minimum is marked with a star while the local 
minima are denoted with a dot. 
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FIG. 2. [CC] and [NC]/[CC] predictions at SNO for the allowed regions in the two-neutrino 
mixing scenarios obtained from the global analysis of solar neutrino data at 90 % and 99 %CL. 
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FIG. 4. [CC] and [NC]/[CC] predictions at SNO as for the SMA region in the four-neutrino 
scenario obtained from the global analysis of solar neutrino data at 90 % (lighter) and 99 %CL 
(darker). The dotted line corresponds to the prediction in the case of no oscillations. 
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FIG. 5. [CC] and [NC]/[CC] predictions at SNO as for the LMA region in the four-neutrino 
scenario obtained from the global analysis of solar neutrino data at 90 % (lighter) and 99 %CL 
(darker). The dotted line corresponds to the prediction in the case of no oscillations. 
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FIG. 6. [CC] and [NC]/[CC] predictions at SNO as for the LOW region in the four-neutrino 
scenario obtained from the global analysis of solar neutrino data at 90 % (lighter) and 99 %CL 
(darker). The dotted line corresponds to the prediction in the case of no oscillations. 
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TABLES 

TABLE I. Best fit points and GOF for the allowed solutions for the global analysis in the 
framework of two-neutrino mixing. 





Active 


Sterile 




SMA 


LMA 


LOW-QVO 


SMA 


Am 2 /eV 2 


5.0 x 1(T 6 


3.7 x 10~ 5 


1.0 x 10~ 7 


3.9 x 10~ 6 


tan 2 9 


0.00061 


0.37 


0.67 


0.00061 


Xmin 


40.8 


33.4 


37.1 


42.3 


Prob (%) 


27 % 


59 % 


42 % 


22 % 



TABLE II. Allowed ranges of c| 3 c| 4 at 90%, 95%, and 99 % CL for the different solutions to 
the solar neutrino problem. 



CL 


SMA 


LMA 


LOW-QVO 


90 


[0.11,0.31] 


[0,0.43] 


[0,0.32] 


95 


[0,1] 


[0,0.52] 


[0,0.44] 


99 


[0,1] 


[0,0.72] 


[0,0.76] 
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